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Abstract. The two dimensional directed sandpile with dissipation is 
transformed into a (1 + 1) dimensional problem with discrete space and 
continuous 'time'. The master equation for the conditional probabil- 
ity that K grains preserve their initial order during an avalanche can 
thereby be solved exactly, and an explicit expression is given for the as- 
ymptotic form of the solution for an infinite as well as for a semi-infinite 
lattice in the horizontal direction. Non-trivial scaling is found in both 
cases. This conditional probability of the sandpile model is shown to 
be equal to a A'-spin correlation function of the Heisenberg XX spin 
chain, and the sandpile problem is also shown to be equivalent to the 
'random-turns' version of vicious walkers. 



1. Introduction 

Non-equilibrium dynamic systems have been for some time of consid- 
erable interest as they can exhibit critical behaviour in close analogy with 
systems at thermal equilibrium. A certain class of such dynamic systems, 
various sandpile models [IjiTj, have become a standard framework when an- 
alyzing self-organized criticality [8l|9], i.e. when the dynamics of the system 
inevitably drives it to a critical state independent of the initial state. De- 
spite the extensive work on these systems, it is only fairly recently that a 
more detailed understanding of problems like when exactly sandpile models 
exhibit self-organized criticality, or what are the possible universality classes 
of their critical behaviours, have begun to emerge. 

Most of the work so far on sandpiles has thus concentrated on properties 
such as e.g. the average duration of avalanches and their size distribution, 
which both exhibit scaling in a critical state. However, there may well 
be for example interesting many-particle correlations in sandpiles, which 
likewise exhibit scaling. If one considers particles with non-intersecting tra- 
jectories, interesting connections with problems like vicious walkers [lOllllj 
would probably arise. Non-intersecting Brownian walkers have also been of 
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very recent interest, and they as well seem to display corresponding scaling 
properties [12]. 

It is the 'fermionic' nature of vicious (and non-intersecting Brownian) 
walkers, which gives rise to scaling of their (asymptotic) survival probability, 
and the related scaling exponent depends in a nontrivial way on the number 
of walkers, as well as on the boundary conditions imposed |13H18j . One 
would thus expect that suitably defined conditional (many-particle) proba- 
bilities of particles with non-intersecting paths in sandpiles should exhibit 
rather similar properties. If the trajectories of vicious walkers and non- 
intersecting Brownian particles are equivalent to 'worldlines' of free fermions, 
one would expect in addition that these sandpile probabilities can be trans- 
formed into a problem of free fermions, or, equivalently, into one of spin 
{S = 1/2) chains. 

In order to address these questions, we consider in this paper the prob- 
ability that K particles of a two dimensional (2D) directed sandpile of ref. 
[1] preserve their initial order during an avalanche. First we reformulate 
the (Abelian) sandpile model such that it becomes one in 1-1-1 dimensions, 
and assume for the sake of generality that it is dissipative, i.e. that the 
number of 'grains' is not conserved in the topplings of unstable sites. We 
also consider two different boundary conditions, an infinite system (in the 
horizontal direction) and a system with an absorbing boundary at the origin 
(a 'semi-infinite' system). As it is well known by now, the model is criti- 
cal only at vanishing dissipation [6l[T9]. We derive an exact analytic form 
for the probability, and show that non-zero dissipation introduces an expo- 
nential cutoff in its asymptotic form that also includes a power law with a 
scaling exponent that depends nonlinear ly on K, and is different for the two 
boundary conditions. 

We also show that this probability is equal to the partition function of 
K vicious walkers, more precisely the 'random-turns' version of such walk- 
ers [101114] . The former probability is thus the generating function for the 
survival probabilities of the walkers. The scaling exponents of the sandpile 
probability are not those of the survival probability of the 'lock-step' vi- 
cious walkers, although the two problems are intimately connected. Finally 
we show that the sandpile probability is equal to a correlation function of 
the Heisenberg XX spin chain. This establishes the relation of both the 
sandpile problem and the walker problem to a problem of free fermions, as 
the Heisenberg XX chain is equivalent to free fermions via a Jordan- Wigner 
transformation. 

2. Abelian sandpile model 

2.1. Discrete model. A 2D directed (Abelian) sandpile model on a 
lattice (see, e.g., [2j) is constructed such that to each site (j, n) an integer 
height variable (number of grains) 2(j,n) is assigned. The site has a threshold 
height z^j below which it is stable. The dynamics of the model consists 
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of two steps. First, we choose a site {j,n) at random and add one grain 
to it, i.e., Z(^j^n) ^ ^(j,n) + 1- For z^^j^^) ^ ^{j n)^ (i;"-) becomes unsta- 
ble and its grains are distributed among the 'downhill' neighbouring sites. 
In the following we will use the notation by which the locations of lattice 
sites in the horizontal direction are labelled by j, k or I, and by n in the 
downhill direction. By n we can equivalently denote the number of steps in 
a cascade of toppling processes. In a toppling at site (j, n) grains are thus 
distributed to sites + + and (j — l,n + l). By supressing the n labels 
(understanding that two adjacent columns in the lattice are connected in a 
toppling and that there is no n dependence) we can express a toppling in 
the form 

^ ^3 - ^ij^ (2-1) 
in which the elements of the toppling matrix A satisfy Ajj > 0, and Aij < 
for I ^ j. The condition Aij ^ for every / guarantees that no grains 
are created in the toppling process. Without loss of generality we can put 
Ajj = z^j . The allowed number of grains in a stable site (j, n) is now 
1,2, . . . , Ajj — 1. The sites {j,n) such that ^j Aij > are called dissi- 
pative. Boundary sites are always dissipative so that grains can leave the 
system through the boundaries. After an initial toppling at a site, neigh- 
bouring sites can also become unstable, and sites are kept on relaxing with 
parallel updating until all sites are stable. In this way an avalanche of top- 
plings is generated. Existence of dissipative sites ensures that all avalanches 
terminate in a finite time. 

Assume now that all lattice sites are initially in a stationary state (i.e. 
are stable): ^(j.n) = ^fjn) ~ If add a grain at a randomly chosen 
site (/,0), and make site {j,n) dissipative such that the system returns to 
a stationary state after the extra grain disappears from this site (i.e. after 
n steps). The conditional probability Gji{n) that an extra grain is at site 
(j, n) satisfies the equation 

Gjiin) = i {G,+iKn - 1) + Gj_u{n - 1)} , (2.2) 

with the initial condition Gji{0) = 6ji. Since we consider only symmetric 
topplings, the conditional probability also satisfies Gji{n) = Gij{n). It is 
easy to verify that eq. (12. 2p is the same as the equation for the corresponding 
probability expressed in the conventional 'light cone' coordinates, eq. (5) in 
ref. [1]. 

2.2. Continuous 'time' model. We can also express eq. (j2.2p in the 
form 

Gj^i{n + 1) - Gj,i{n) = ^ {G,+i,,(n) + G,_i,,(n) - 2G,y(n)} . (2.3) 

Consider now a process in which the discrete number of steps is replaced 
by a continuous parameter that will be called 'time' in the following. Let 
Pji{t) be the conditional probability that a grain is at a horizontal location 
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j at time t after an arbitrary number of steps since it was dropped at a 
horizontal location / at t = 0. Transforming eq. (|2.3p into such a continuous 
time we obtain that, during a short time interval dt, the probability Pji{t) 
changes such that 

Pj,i{t + dt) - Pj,i{t) = ^ {Pj+i,iit) + Pj-i,iit) - 2Pj^iit)} dt, (2.4) 
which leads to the master equation 

^^P^i{t) = -^Y.A,kPki{t) (2.5) 

k 

with the toppling matrix 

Ajk = 26jk - {Sj+i,k + Sj-i^k)- (2.6) 

This toppling matrix means that, as above, at each toppling two grains 
are removed from the site and distributed to its nearest-neighbour downhill 
sites. We consider here only symmetric topplings, Ajk = A^j, and thus 
Pjk{t) = Pkj{t). The initial conditions are Pjk{0) = 6jk- For a model of 
N sites in the horizontal direction, the lateral boundary elements of the 
toppling matrix can be defined such that Ao,i = An,n+i = 0, and hence 
the boundary sites j = 1 and j = N are always dissipative as required. 

Notice that the continuous time is not a simple continuum formed by the 
discrete variables n, but Pji{t) includes processes with all possible numbers 
of steps. In fact function Pji{t) can be considered as the generating function 
of the conditional probabilities Gji{n) as we find that 

4.n 

e'p^,{t) = Y,G,i{n)-. (2.7) 

n=0 

The expected number of topplings at site j in an avalanche resulting from 
a perturbation (adding a grain) at site / is given by Tji{^) = X^^o^iK'^)- 
A Laplace transform of the conditional probability 

oo 

r,K/) = j e-'fPji{t)dt, (2.8) 



is the Green's function of the master equation eq. (j2.5p . It is easy to verify 
that Tji{0) satisfies [2j the condition ^jk^kii^) = ^ji- 

The master equation eq. (j2.5p can easily be solved for the toppling 
matrix of eq. (j2.6p with the initial condition Pji{0) = dji- Consider first the 
case of an infinite lattice in the horizontal direction such that — c« < j,l < 
oo. In this case we find that 

Pj-^'^\t)=e-'li_,it), (2.9) 

where Ij{x) is a modified Bessel function. Asymptotically, for large t, the 
conditional probability for a single grain thus behaves as 

pj-oo,oo)^^^ ^^_l/2^ (2.10) 
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There is a pure power law so that the duration of avalanches scales with an 
exponent C(_oo,oo) = 1/2- This exponent coincides with the known result for 
2D directed sandpiles [T| as it should. 

We can analyze the effect of boundary conditions by introducing an 
absorbing boundary at the origin. To this end we first recall the solution for 
a finite lattice of N sites (see e.g. |19) ) for which the boundary conditions 
are Pji (t) = for j,l = and j,l = N + 1. In this case one has 

P,(t) = e~*^^-sin-^sin-^, (2.11) 

J'^ ^ 7V + 1^ N + l iV + 1' ^ ' 

k=i 

where the spectrum is of the Bloch form, 

Ek = l- cos . (2.12) 

iV + 1 ^ ^ 

In the limit N ^ oo the sum in eq. (j2.1ip can be replaced by an integral, 
with the result 

TT 

pj°'°°)(t) = 1 J e-*(i-™«^) sin(/x) sin(jx) 


= e-*[/i_,(t)-/i+,(t)]. (2.13) 

The asymptotic behaviour for large t of the conditional probability is now 
given by 

Pj°'°")(t) ocr3/2. (2.14) 

The scaling exponent indeed depends on having a boundary at a finite dis- 
tance: C(o,oo) = I = C(-oo,oo) + 1- As expected, the same exponent has been 
found for the scaling of avalanche sizes with the corresponding boundary 
conditions [201I2T]. 



3. Multiple-grain correlations 

3.1. Master equation. Having established that our master equation 
method indeed reproduces previously known results, we turn now to a more 
interesting problem of correlations between multiple grains during 'avalanche 
dynamics'. To this end, let us address the following problem. Consider the 
same lattice as above with all its sites in a stationary state: -Z(j,n) = n)~^^ 
and add K grains at randomly chosen K horizontal locations: li > I2 > 
■ ■ ■ > Ik- The toppling rules are the same as above, at each toppling two 
grains are removed from the toppling site j, and a grain can jump to each 
of the two nearest-neighbour sites in the downhill direction. However, if 
■^0'.") ~ ^(j±i,n+i) = 0) site {j, n) cannot topple. The probability that the 
additional grains will be at dissipative sites ji > 32 > ■ ■ ■ > Jk at time t 
(after an arbitrary number of topplings) satisfies a generalized version of eq. 
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(j2.5p . namely 



d I ^ 

^-^ilv,iA';'lv,'A' (*) ~ 2 [Pj-i,--;jK;h,---,lr-l,lr + l,lr + -L,---,lKi^) 

r=l 

+ Pjl,...,jK,h,-,lr~l,lr-l,lr+l,-,lKii^)] " ^ ^jl ,■ ■ ■ ,j K A ,■ ■ -,1 K i^) ^ (^-l) 

supplemented by the condition i-*;i,...j^;Zi,...,Zj^(t) = 0, if jr = jr+i, for all 
r = 1, . . . , K — 1. The solution to this equation is given by 

Pn-jKA-iAt) = ^ det {P,vJt)}, (3-2) 

where Pji (t) is the one-grain conditional probability which satisfies eq. ()2.5p 
with the same boundary conditions as the solution of eq. (13. Ih . 

As for the single grain, in the multi-grain case the continuous conditional 
probabilities Pji....,jK;h----,iKi^) generating functions of the discrete ones, 
^iivjV;'i,---,'A('^)' '^^ ^^'^ ^^^^ 

e-^^Ph,-,3K;h,-4At) = '^C!j^^,„jj^.i^^,„^ij^{n)——. (3.3) 

n=o 

The discrete probabilities satisfy the equation 

K 



1 V— -\ 

=1 

+ Gji,...j^_,jv_ij^+i,...j^;;,,...,z^(n - 1)}, (3.4) 



r=l 



supplemented by the condition Gj^,,,,^jj^-i^,,,,^ij^{n) = 0, if jV = jV+i, for all 



3.2. Infinite lattice. Let us now consider the asymptotic behaviour for 
i — >■ oo of the above multi-grain conditional probability. We first consider 
the case of an infinite lattice in the horizontal direction when the one-grain 
probability is given by eq. ()2.9p . Using the integral representation for the 
modified Bessel function, we arrive at the expression 

TT TT 

^■""rl , {t) = t\k fdxi - - [ dxK e-*^^=i(i-'=°^^'") 

— TT — TT 

X det /e*('»->)^'-l . (3.5) 

Making use of the symmetry of the integrand with respect to permutations 
of integration variables xi, . . . ,xk, the determinant in this expression can 
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be transformed such that 



K 



det (e*(^^->)^'-| ^ det I e^^^^''' | TT ' 

r=l 

det je"*^'^^'-) det \e'^^^A. (3.6) 

The two determinants above can be represented in terms of Schur functions 
(for a survey on Schur functions see e.g. [22j ) : 



sxixi,X2, . . .,xk) ■■ 



s 



det (x^+^-^•) n i^s-x,)-\ (3.7) 



l^s,k^K 

l^s<k^K 



where A = (Ai, A2, . . . , Ai^) is a partition of a non-increasing series of the 
non- negative integers Ai ^ A2 ^ . . . ^ Xk ^0. If we consider the case 
jr ^ —K and ^ —K, we find that 



7r TT 



/ (27r)'^i^! J J 



I ixi iX2 ixK\ I —ixi —1x2 

X s\[e ^...,e ^)s^(e ',...,e 



X 



where \r = jr — K + r and = Ir — K + r. 

As t — >• cx) (and js—lr ^ for ah r, s = 1, . . . , i^), the main contributions 
to the above integrals come from near the origin of the integration variables, 
and in leading order we find that 

p(-oo,oo) . s ^ sx{l,l,■■■A)s^^{l,l,■■■A) 
ji,-,jK;h,-,iK^ ) (27r)^/^! 

oo oo 

/ dxi--- / da;^ e"5E^=i^^ JJ (^Xr-Xs)'^. (3.9) 



OO 

X 

-oo 



The prefactor of the integral can be computed (see e.g. [22]) using the well 
known result 

sa(1,1,... ,1) = — — ; , (3.10) 

Ilm=l 



while the integral is the Mehta integral of the gaussian unitary ensemble of 
random matrices [23], which can be explicitly evaluated: 

-oo -co l<r<s^K 

(3.11) 

We thus find that, as t — )• oo, in leading order the multi- grain conditional 
probability is given by 



with the scaling exponent 



7 = ^ (3.13) 



and the amplitude 

nisgs<r^_ft:(^'' ~ ^s){jr — is) 



Aj^,...,jK\h,-,lK - K-1 '. • > 



3.3. Semi-infinite lattice. Let us consider the conditional probability 
in the presence of an absorbing boundary at the origin. As in the one-grain 
case, let us start with a finite lattice of N sites in the horizontal direction. 
Substituting eq. (|2.1ip into eq. (|3.2p . we find that 



Pji,-,jK;h,-,lKi'^) (]\f_i_-[\K X/ 



K N ^ _ K 

^ ' fci = l kK = l 



X det I sin sin I , (3.15) 

where is given by eq. (|2.12|) . The multi-grain conditional probability for 
the semi-infinite lattice follows from this result by taking the large N limit; 
the resulting expression is similar to eq. ()3.5p . but with a determinant that 
now contains sine functions instead of exponential functions: 

7r TX 

p(0'~). (t) = ^ [ dxi--- [ dxK e-*S^=i(i-'=°^^-) 

Jl,--;jK;h,--;lK^ ' J J 

— TT — TT 

X det |sin(7Vxr) sin(/sa;r)l • (3.16) 

Again, using the symmetry with respect to permutations of the integration 
variables xi, . . . ,xk, we can transform the determinant in this expression 



such that 



K 



det {sin(jVa^r) sm(/sX.r)} — > det {sm(/sXr)} IT sin(j' 



r=l 

(3.17) 

Using the character of the irreducible representation corresponding to a 
partition A of the symplectic Lie algebra, 

Spx{Xi,X2,...,XK) - — K_k+^ -(K-k+1). ' (3.18) 



we can express eq. (j3.16p in the form 



X det {sinsi,} spA(e"',e™,...,e"'<) 

X sp^(e^^\e"\...,e^^'^), (3.19) 

where A,- =jr — K + r — 1 and /Ur = /r — -^ + '' — 1- The determinant in the 
above integrand can be evaluated using the identity (for a proof, see 



K X- -X X- +X 

det {sinsxr-l = 2^*-^^^'' TT sinxr- TT sin sin^^ . 

(3.20) 

We finally obtain for the conditional probability the expression 



p{o,oo) {t) = - dx^■■■ / dxA- e-*S-=i(^-^°''^'") 



22K(K-1) 

—TV — TT 

n. 2 TT ■ 2 -^j ■ 2 -^i "T ''"fc 

sm Xr YY sm ^ sm ^ 



X 



spA(e*"^ , e'"' , . . . , e^"^)sp^(e^^'i , e^"^ , . . . , g^-^c). (3.21) 



In the limit t — )• oo we can approximate the integrals in the above 
expression with the integrals 



oo oo 



nLi(2"^ 



\^3<kii,K j=l 

- (27r)^/2t^{2i^+i)72- (^-^^^ 

For a proof of eq. ()3.22p , see [23] . We find thereby for the leading asymptotic 
term of the generating function 

Here the scaling exponent is given by 

i^(2i^ + 1) 
^ = 2 

and the amplitude is 



(3.24) 



nLi(2m) 



An,-,jK^iu-,iK = ^7^^t-sPa(1, . . . , l)sp,{l, . . . , 1), (3.25) 
in which 

.PA(1,...,1)= n ijr-fs)U ^'%T^^f - (3-26) 

l^r<s^/^ m=l ^ ^ 

A similar expression can be found for sp^{l, . . . , 1) with the jV's replaced by 

We thus find that the scaling exponent of the multi-grain sandpile prob- 
lem considered here is not equal to the one found previously for the 'lock- 
step' version of vicious walkers, for which 7 = K{K — l)/4 [10','11]. Instead, 
our sandpile problem corresponds to the 'random-turns' version of vicious 
walkers |10tll4j. The connection to the 'random-turns' version of vicious 
walkers is discussed in more detail below. 



4. Connection to vicious walkers 

4.1. Heisenberg chain. Before addressing the relation between the 
above sandpile problem and the 'random-turns' vicious walkers, we first 
outline its relation to free fermions. It turns out that it has a straightforward 
connection to the Heisenberg XX spin chain that can be mapped, as is well 
known, to a free fermion problem by the Jordan- Wigner transformation. 

The Heiseberg XX chain has the Hamiltonian 

i,k 
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where summation is over all lattice sites, and 

Ajfc = 5i^k+i + ^i,k~-i- (4.2) 

We use the standard notations crf,af for Pauli spin operators that satisfy 
the commutation relations 

[a+,a,] = a^6,k, , a±] = ±2at6,,, (4.3) 

and have in addition the properties 

{atf = 0, (af )2 = 1. (4.4) 

In what follows we use the fact that the ferromagnetic state with all spins 
up, I lT ) = ®i 1 1 )i) which satisfies o"^ | ft ) = for all k and normalized such 
that ( ft I lt ) = 1, is annihilated by the Hamiltonian, 

H\^) = 0. (4.5) 

Our aim is to study the 'temporal' evolution of states with a finite number 
of down spins, which can be constructed by acting with operators a~ on the 
state I ft ) . We thus consider the matrix elements 

Fj,,...,j,;i„...,i,{t) = I a+ • • • CT+ e-*^a- • . . af^ (4.6) 

Parameter t will play the role of 'time' in the context of the sandpile model. 

Before proceeding with the general case, let us first consider the case 
K = 1, i.e. the 'temporal' evolution of a single reversed spin. Differentiating 
the function Fji{t) = ( ft | cr^e^*^cr^^ |ft) with respect to t and using the 
commutation relation 

[^'^1 = -IE^^>^^H = -r'A^U + '^t+i) (4-7) 

k 

together with the property ( ft | (t| = ( ft | , we find that 

j^F^iit) = -{it\alHe-'^a-\t) = ^( ^ I + <+i)e-* V I ^ )• (4-8) 
Hence the correlation function eq. (j4.6p for K = 1 satisfies the equation 

j^Fjiit) = \ iF,+u{t) + F,-iKt)) . (4.9) 

Similarly, by commuting H with , a difference equation similar to eq. 
(j4.9p follows, but for subscript j with fixed subscript /. Both equations are 
subject to the initial condition Fji(0) = 5ji, and to boundary conditions that 
depend on the type of the lattice: for the semi-infinite lattice Fji = for 
j, I = 0, while for a finite lattice Fji = for j,l = and j,l = N + 1. 

As sl result, comparing eqs. (14. 9p and (12. Sh together with their initial 
and boundary conditions, we find that the one-spin correlation function of 
the Heisenberg XX chain is equal, modulo a trivial factor, to the one-grain 
conditional probability of the sandpile model, i.e. Pji{t) = e~^Fji{t). 
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Let us now consider the case of general K. Differentiating eq. (|4.6p with 
respect to t, taking into account the differential property of the commutation 
relations, 

K 

[a+(7+---(7+ F] = \(7t,H\(yt ■■■'yt , (4.10) 

k=l 

and applying the commutation relation eq. (|4.7|) . we find that 
d 1 ^ 



+ ^ii,...j,_ijV-ljr+i,...jK;h,...,«A-(Oj- (4-11) 

A similar equation can be found with respect to subscripts Ir with the jV's 
kept fixed. The initial condition is -?ji,...jk;'iv,'a'('-') = '^ji^i ' ' ' ^okIk- '^^^ 
correlation function also satisfies the conditions Fjx,..;jK\h-- -,iKi^) = if 
Ir = Is or = (r, s = 1, . . . , K), which follow from the nilpotency of the 
Pauli spin operators, eq. (j4.4p . 

It is evident that the differential equation eq. (|4.1ip which the correla- 
tion function eq. (j4.6p satisfies, coincides with eq. (j3.ip for the multi- grain 
probability up to a trivial 'diagonal' term. It is also easy to verify that the 
solution of eq. (j4.1ip can be expressed in a determinant form, 

F,,,...,,^;i,,...,iAi) = , det , (4.12) 

where Fji{t) are the one-particle correlation functions satisfying eq. (j4.9[) . 
We thus find that 

Pji,...,jK;h,-,lKi*) = (i~^^Fh,...,jK;h,-,lKi^)- (4-13) 

Hence, the multi-grain probabilities in our sandpile problem are also multi- 
spin correlation functions in the Heisenberg XX spin chain. 

4.2. Quantum trajectories. We are now in the position to establish 
an explicit relation of our sandpile model with random walks. To this end 
we exploit the known connection of the Heisenberg XX spin chain with the 
random-turns vicious walk model [25ti27j . Our starting point, which follows 
from eqs. (|4.13p and (|3.3p . is that the discrete multi-grain probability can 
be expressed as a matrix element: 

By evaluating this matrix element, it can be shown that the discrete proba- 
bility of the sandpile model is equal, modulo a simple factor, to the number 
of paths of random-turns vicious walkers subject to the same boundary con- 
ditions. 
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We start again with the single-grain case in which the discrete probability 
Gji{n) is the one-spin matrix element ( ff | (T^(— -ff)"(T^^ | f|~ ). Taking into 
account eq. (|4.5p . we can write 



fci 

and repeating the procedure for the nth power, we find that 

^Vl^) = ^ E ^kr.k„.,---^k,k,Ak,ia,Jil). (4.16) 

Multiplying this expression from the left by ( j and using the orthogo- 
nality of the spin states, ( f|~ [ (^j^ctj^ \ = ^ju we find that 

(^|a+(-i7)VI^) = ^ A,k„^,---Ak,kAk,i- (4.17) 

h\ , . . . ,/?n — 1 

The sum in eq. ()4.17p can be interpreted as one over all possible quantum 
trajectories (lattice paths) of n time steps of a particle (corresponding to 
the down spin state) from site / to site j subject to the boundary conditions. 
In this interpretation A appears as the transfer matrix. From eq. (|4.2p it 
follows that we deal with a random walk on a lattice. Denoting by Vnil j) 
the number of all admissible paths of n steps from site I to site j, we have 

Gji{n) = ^Vnil^j). (4.18) 

Let us now consider the general case. Acting with the Hamiltonian H 
on the state c^^cr^^ • • • crf^ lit), for which we assume that li > I2 > • • • > Ir, 
we find that 

K 



r=l m 

~ ~2 ^mi,...,m/f ;Zi,...,iK'-''mi ■ ■ ■ '-''mji I 1T ); (4.19) 

mi,...,mii 

where 

K 

'^mi,...,mjf;li,...,lj( ^ ^ ^mih ' ' ' ^rrir-ilr-i-^'mrlr^Tnr+ilr+l ' ' ' ^it^kIk (4.20) 



r=l 



for mi > m2 > • • • > rriK, and r„^,...,m^;/i,...,z^ = for = m^+i (r = 
1, . . . , -ftT — 1). Interpreting T as a transfer matrix, we find for the multi-spin 
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matrix element an expression in terms of a matrix element of the nth power 
of this transfer matrix, 

I < ■ ■ ■ ^tJ-^r^h ■ ■ ■ I ^) = ^iTln,-M,...,iK, (4.21) 

which generalizes eq. (I4.16P to the case of K random walkers. 

The presence of just single factor A in each term of the transfer matrix 
eq. (j4.20p implies that, at each time step, only a single walker moves out 
of the total K. Thus the above matrix element of the nth power of T gives 
the number of all lattice paths of n steps made by K random-turns vicious 
walkers. We recall that, in the random-turns vicious walkers model, at 
each time step only a single randomly chosen walker moves one step to the 
'left' or one step to the 'right' with the constraint that two walkers cannot 
occupy the same site. These random-turns vicious walkers are different from 
the more common lock-step vicious walkers all of which, at each time step, 
must move left or right with the same constraint that two walkers cannot 
occupy the same site [lOj. 

Denoting by Vnih, ■ ■ ■ Jk ^ ji, ■ ■ ■ jJk) the number of all admissible 
configurations in which the K walkers are initially located on the lattice 
sites li > I2 > ■ ■ ■ > Ik, and have after n steps arrived at the positions 
ji > 02 > • • • > Jk, we find the result 

Gn,-,jK;iu-,iKir^) = J^'^nih, ...Jk^Ji,.. -Jk), (4.22) 
where RT stands for random-turns vicious walks. 



4.3. Large n limit. Having established a connection of our sandpile 
model with the random-turns vicious walkers, it is natural to consider the 
large n limit of discrete conditional probabilities. This is useful for a direct 
comparison with the random-turns walkers, for both an infinite and a semi- 
infinite lattice in the horizontal direction. 

For definiteness we consider here an infinite lattice in the horizontal 
direction; a semi-infinite lattice can be considered similarly, and below we 
outline the results for both cases. Using the relation between the continuous 
and discrete conditional probabilities, see eq. ()3.3p . we find from eq. ()3.5p 
in the case of an infinite lattice the representation 



^(-00,00) , s 



(2vr) 



-YJ^ / rfxi • • • / dxK ( ^ COSXjn j 

— TV —TV ^ 

X det (e*('--^'-)^' | . (4.23) 



We are interested here in the large n limit with the /^'s and jV's kept 
fixed. In order to apply the standard saddle-point approximation, we ex- 
press the first factor of the integrand in the above equation in the form 
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exp {n log {J2m '^^^ ^m)}; and obtain thereby the fohowing system of saddle- 
point equations: 

^k"^' =0, r = l,...,K. (4.24) 

It is evident that the solutions to this system of equations satisfy sinx,- = 
0(r = 1, . . . with the restriction that 'Ylm 

COS XfYi ^ 0. Requiring that 

the matrix of second derivatives 

cosxr- ^ sinxrsinxs r.r-\ 
brs — ; -o (4-^5) 




is a negatively definite matrix for the solution of eq. ()4.24p . we find that 
the steepest descent corresponds to the solution for which cosx,- = 1 (r = 
1, . . . , K) i.e. the main contribution to the integrals in eq. ()4.23p comes from 
near the points = (r = 1, . . . ,K), similarly to the case of continuous 
conditional probability considered in sect. 3.2. 

Therefore, replacing the first factor of the integrand in eq. ()4.23p by its 
approximation near the origin of the integration variables, i.e. 

f'j^cosx^'j ai^"exp|-^ J^x^j, (4.26) 

\m=l / I. m=l ) 

and transforming the second factor of the integrand as in sect. 3.2, we find 
that as n — >■ oo the leading order form of the discrete multi-grain probability 
in the case of an infinite lattice can be expressed as 

( — 00,00) / N 



^(,-00,00; 

00 00 

X f dxi--- f dx/^ e-("/2^)^™=i^™ Yl i^r-xsf. (4.27) 

-L -00 li:r<si:K 

Clearly this expression is the same as eq. (13. Oh with t replaced by n/K. 
It is easy to check that a similar result is valid for a semi-infinite lattice, 
now using the procedure of sect. 3.3. Hence, the leading terms of the large 
n limits of discrete multi-grain probabilities can be obtained from those of 
the large t limits of the continuous probabilities simply by the replacement 
t H> n/K. 

We have thus shown that, as n — )• 00 for fixed Ir's and jV's, the discrete 
multi-grain conditional probabilities scale as 

Gn,...dKM,-,iKin) ~ Bj^,-,jK;h,-,iK ^"^ (4-28) 

with 

Bji,...,jK;h,...,iji = K'^Aj^^^jj^-i^^^^ij^, (4.29) 

where the exponent 7 and amplitude given by eqs. (j3.13p 

and (13.14p . respectively, in the case of an infinite lattice, and by eqs. (|3.24p 
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and (|c{.25|) . respectively, in the case of a semi- infinite lattice. The discrete 
forms of the conditional probabilities scale exactly as the continuous ones, 
as they should. These results, in view of eq. (j4.22p . are also in agreement 
with the known scaling properties of the random-turns vicious walkers [13]. 

5. Conclusion 

In conclusion, we showed that the probability of multiple grains in a two 
dimensional directed sandpile to preserve their order during an avalanche, 
displays non-trivial scaling properties that are similar to those of the 'ran- 
dom-turns' version of vicious walkers. The conditional probability in the case 
when the downhill direction of the lattice was transformed into a continuous 
variable ('time'), which we found as the solution to a master equation, was 
found to be the generating function of directed lattice paths in the origi- 
nal lattice. This continuous time conditional probability then provided the 
connection with the Heisenberg XX spin chain as it was found to be the 
same as the corresponding many-spin correlation function of the Heisenberg 
chain. Notice that connection with spin- 1/2 variables was based on having 
a directional lattice with the chosen initial condition. After a toppling on 
'row' n there were only sites with a critical or subcritical (critical minus 
one) number of grains on row n + 1, i.e. there were effectively only two 
possible states per site. For other situations connection would possibly be 
to spin systems of higher spin. Connection with the Heisenberg chain also 
established a relation between sandpile models and free fermions. Similarly 
to the sandpile problem, the spin correlation function was shown to be the 
generating function of 'vicious' quantum trajectories of spins. 
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